Electric field effect on superconductivity at complex oxide interfaces 
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We examine the enhancement of the interfacial superconductivity between LaAlC>3 and SrTiOs 
by an effective electric field. Through the breaking of inversion symmetry at the interface, we show 
that a term coupling the superfluid density and an electric field can augment the superconductivity 
transition temperature. Microscopically, we show that an electric field can also produce changes 
in the carrier density by relating the measured capacitance to the density of states. Through the 
electron-phonon induced interaction in bulk SrTiOa, we estimate the transition temperature. 



I. INTRODUCTION 

Superconductivity is typically thought to be produced 
by the pairing of electrons below a critical temperature 
T c »i The overall pairing interaction is dependent on a 
subtle balance of the repulsive Coulomb interaction and 
any attractive interaction supplied by the exchange of 
collective excitations in the systemjL 2 - That balance de- 
pends, in a complex way, on the number of carriers within 
the materiali 2 - - — This effect becomes all the more impor- 
tant in the case of spatial inhomogeneity, e.g. at an inter- 
face, where the spatial distribution of the carrier density 
will be important Therefore, it is of great interest to 
have a tuning parameter like an external electric field 
producing a change to the carrier distribution and effec- 
tively maximizing the transition temperature T c . 

In recent years, oxide hctcrostructurcs have attracted 
much attention due to a rich spectrum of unexpected be- 
havior discovered at interfaces^ Specifically, the emer- 
gence of a 2-D electron gas, that becomes superconduct- 
ing at around 200 mK, at the epitaxially grown interface 
of the insulator SrTiC>3 (STO) and a few atomic lay- 
ers of LaA103 (LAO) (illustrated in Fig. [IJ has sparked 
an intense search for the mechanism of superconductiv- 
ity as well as for improved materials <2r— Studies investi- 
gating the electron mobility have demonstrated that in- 
trinsic doping may be key to understanding the under- 
lying mechanisms^ While the superconductivity occurs 
at much lower temperatures than in high-T c supercon- 
ductors, this finding shows the importance of investigat- 
ing interfacial physics in heterostructures. Even though 
there are limitations due to lattice matching and strain 
within the materials , 13 ! 14 the coupling between these ma- 
terials may provide information that is critical to the un- 
derstanding of these heterostructures. 

Further studies of the LAO/ STO system have shown 
that T c can be increased by an applied electric nek h 15 i 16 
Bulk STO has been shown to become superconduct- 
ing with doping by oxygen vacancies and by niobium 
(Nb) at about the same temperature as observed in the 
LAO/STO systemii Furthermore, STO has been found 
to exhibit interfacial superconductivity when brought in 
connect with an oxide gel in the presence of an external 



electric fieldJ^ These studies show that T c can be in- 
duced and/or enhanced by the creation of carriers pulled 
from a source (i.e. LAO, oxide gel, or doping) i 15 i 16 ' 18 

The nature of the superconducting state at the in- 
terface is a subject of ongoing discussion. One pos- 
sibility is that the superconducting state is conven- 
tional electron-phonon driven, but there is also a grow- 
ing list of possible unconventional states including 
magnetism driven onesi 19 ' 20 Indeed recent data from 
scanned probes point to a complicated spatial pattern 
of coexisting superconductivity and magnetism at the 
interfaces^ Early reports on ferromagnetic order coex- 
isting with superconductivity^ 2 - have been independently 
confirmed^ For an assessment of the present situation 
see Ref. [24 1. 

For the following, it will be important to know the 
thickness of the mobile electrons layer at the interface. 
Early estimates of the 2-D density of mobile carriers 
have yielded values much below the theoretical predic- 
tion of 0.5 electrons per unit cellj 2 ^ Apparently, most of 
the added charge is not mobile and is trapped within a 
few unit cells of the surface, as observed by high-energy 
photoemission experiments^ while the mobile electrons 
may cover a region of thickness 10 nm or more. 

The recent experiments 2 ^ suggest that the spatial sepa- 
ration of superconducting and ferromagnetic regions may 
allow spin singlet pairing to coexist with ferromagnetic 
domains. The spatial extension of the density of mobile 
carriers normal to the interface over several tens or even 
hundreds of unit cells speaks against a prominent role of 
the on-site Coulomb interaction and therefore a pairing 
mechanism induced by magnetic fluctuations. All of this 
points to a conventional s-wave superconducting state. 

In this study, we investigate i) the effect of an elec- 
tric field applied externally on the 2-D and 3-D carrier 
density. We calculate the variation of the 2-D carrier den- 
sity, capacitance, and the 3-D carrier density depth pro- 
file with increasing electric potential for LAO/STO het- 
erostructures. From changes of the overall carrier density 
at the interface, we deduce that the nonlinear capacitance 
of the layer at the interface implies an energy depen- 
dent Density of States (DOS); ii) using modifications of 
charge carrier density and DOS as a tunable parameter, 
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FIG. 1: An illustration of LaAlOs/SrTiOs interface, which 
details the interface region. The interfacial (lighter color) area 
is a small region from — Sz to Sz. The bulk (darker color) re- 
gion is that volume that is essentially unaffected by the inter- 
face interactions. For the purpose of this study, we examine 
the effect of an applied electric field produced throughout the 
sample. Therefore, all calculations are to mimic this general 
setup. 



we calculate changes in superconducting transition tem- 
perature T c . We assume the pairing to be conventional 
s-wave pairing as identified in detailed studies of doped 
bulk STO ii£ Using this model assumption, we show that 
the T c for the band insulator STO can be modified by 
changes in charge density, and we find that T c does ex- 
hibit a maximum at a certain carrier density determined 
by the peak in the capacitance. We show a plausible con- 
nection between the increase in carrier density and T c due 
to an applied electric field. This observation could lead to 
a possible mechanism for a tunable functional response: 
dramatic changes in the superconducting properties with 
the application of external fields. 



II. FIELD EFFECT ON THE 
GINZBURG-LANDAU FREE ENERGY 
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FIG. 2: Model calculation for the shifting electrostatic dif- 
ference potential (EDP) due to an applied electric potentials 
ranging from V to 10 V using a semi-empirical, extended 
Hiickel model for the LAO/STO interface. The model 
shows the diffusion of the potential (approximately I unit 
cell per 2V) from the LAO region into the STO region (the 
arrows are to help guide the eye). In the simulation, the last 
two unit cells are used as the electrodes (LAO - cathode and 
STO - anode) producing an overall electric potential over the 
sample region. 



The breaking of inversion symmetry at the interface 
of heterostructures introduces the interface normal vec- 
tor as a preferred direction of the system.— The normal 
component of the electric field F z is therefore compat- 
ible with the symmetries of the free energy, allowing a 
coupling term of electric field and superconducting den- 
sity, F z \ip\ 2 , with the superconducting order parameter 
ip. This is a simple consequence of screening effect on the 
charge density within the superconductor^ The charge 
distribution at the interface would be dependent on the 
distance from interface, typically on the order of nm. On 
the other hand the relevant superconducting length scale 
is set by a fairly large coherence length £/2, typically on 
the scale of tens to hundreds of nm4 Therefore, we may 
use a Ginzburg-Landau (GL) functional in our approach 
when we address the coupling between SC and electric 
degrees of freedom^ 

Since the interface has the additional constraints of 



boundary conditions and strain (as shown in Fig. []}, the 
standard heterostructurc free energy can be broken into 
two main components consisting of bulk (B) and interfa- 
cial (/) regions. The effect of the interface decays off as 
one moves into the bulk. In the case of the LAO/STO 
interface, the total free energy FE is then given by 

/ — Sz p0 
F(z) B sTodz+ / F(z)jstocIz 
-L J-Sz 

(1) 

+ / F{z)jhAodz + / F(z)gi,Aodz, 

JO JSz 

where the depth-dependent free energy for the bulk and 
interface are the given by 

FE{z) B = a B \ij{z)B\ 2 + ^{z)Bt (2) 
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FIG. 3: Simulated capacitance and 2-D carrier density n^d as 
a function of electric potential for the LAO/STO interface in 
Ref. It should be noted that the change in curvature 

denotes an effective change in electron or hole doping. 



and 



FE(z)j = aj^izhl 2 + ^(z),| 4 + g\ V ^)/| 2 

(3) 

+AF z |V(z)/| 2 +7^lV'W/| 2 , 

where asj and /3b, i are the standard GL energy parame- 
ters. In the standard GL formalism, the parameter olb ~ 
as(T — T®), with a > 0, determines the normal (as > 0) 
and superconducting states (as < 0), while /3b is always 
positive and helps stabilize the ground state. 

For the interface region, the gradient and strain 77 
terms are produced by the interface and are typically 
positive, where the gradient produces the boundary con- 
dition for the interface and can be assumed to be in the 
small wavevector limit. For the field interaction, we use 
F = (0,0, F z ). The coefficient A will in general depend 
on the electric field strength. Therefore, the local carrier 
density can increase or decrease depending on positive or 
negative F z , respectively. It should be noted that while 
the electric field is global, only the superconductivity at 
the interface will increase due to the accumulation of car- 
riers from the carrier source. This may be different for 
a high temperature superconductor, but given that STO 
is an insulator, the distribution of carriers is lessened. 

Focusing on the effects of the additional quadratic 
terms, it can be easily shown that the standard GL def- 
initions are changed. The inclusion of the electric-field 
interaction and strain terms induces a shift in T c . Con- 
centrating on the ip 2 terms, we re-write the free energy 



as 

FE(z)f = a/|V(*)i| 2 + AF z |-0(z)/| 2 + 7 ?7 2 |^)/| 2 
-a|^)/| 2 - 

(4) 

where a = otj + XF Z + jr/ 2 = a(T — T c °) and describes the 
superconducting state, when a < 2 -— . Therefore, the 
shifted critical temperature is given by 

T c = T° - ( XcFz \ iW)c ) (5) 

where A c and (jrj 2 ) c are taken at T = T . From this we 
conclude that the electric field has to compete with the 
lattice strain. Even a large lattice strain may be over- 
compensated by a reasonably strong electric field, which 
therefore can still effect the superconducting order pa- 
rameter. However, if the lattice matching is close, then 
the strain can be assumed to have a negligible effect on 
the superconducting state. Recent measurements indi- 
cate the various degrees of lattice mismatch may even 
produce an increase in the superconducting T c ^ 

Thus, the addition of extra quadratic terms to the 
free energy modifies the critical temperature. While the 
strain terms usually produce a negative shift of T c , a suit- 
ably directed electric field could gives rise to a positive 
shift. This may be interpreted as being due to an increase 
in the carrier density of the superconductor. While this 
is shown here only in the phenomcnological GL model, 
we will now turn to a modeling of the charge density in- 
duced by an applied electric field and show how that is 
correlated with T c . 

III. SEMI-EMPRICAL MODELING OF THE 
ELECTROSTATIC DIFFERENCE POTENTIAL 

To examine the change in the electron density at an 
LAO/STO interface, we calculate the electrostatic dif- 
ference potential (EDP) <ff using a semi-empirical, ex- 
tended Hiickcl model. The EDP is described as the dif- 
ference between the electrostatic potential cf> of the self- 
consistent valence charge density and the electrostatic 
potential from atomic valence densities! 30 ' 31 Therefore, 
the EDP is determined by solving the Poisson equation 
V 2 0*=— p*/e r , where is p* is the electron difference den- 
sity. 

Figure [5] shows the calculated EDP for various applied 
electric potentials. This calculation is produced by build- 
ing a three dimensional interface of LAO and STO that 
is 3 x 3 x 8 unit cells for each side. The electric potential 
is produced by simulating a device where positive and 
negative electrodes are place on the ends of the interface. 
For this simulation, the LAO electrode is negative and 
the STO electrode is positive. As the electric potential 
is increased from 0V to 10V, the EDP is observed to dif- 
fuse from LAO into STO. This details the shifting of the 
electron charge density through the interface. To gain 
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a more detailed view of the change in the charge and 
carrier densities, we examine the general electrostatics of 
the interface and compare to known experimental data. 



gularity and is discontinuous at F = —A/B. The singu- 
larity can be removed by expanding out the denominator 
to a quadratic term, producing a continuous function for 
all F, 



IV. CARRIER DENSITY AND CAPACITANCE 

The basic connection of the charge density and an ap- 
plied electric field is given by integrating Gauss's law for 
a dielectric medium, keeping the nonlinear dependence 
of the dielectric function e(F) = eoe r {F) on electric field. 
This function may be determined through the measured 
capacitance at the LAO/STO interface. The change in 
the 2-D carrier density n 2 d in a semiconducting channel, 
induced by an electric field is then given by 



n 2d (F) = ^- l F e r (F)dF 
e Jo 



(6) 



where F is the effective electric field within the sample, 
eo is the permittivity of free space, t r (F) is the electric 
field dependent relative permittivity, and e is the basic 
unit charge. The field dependence of the relative per- 
mittivity for STO has been shown to be e r (F) STO = 
1/A (l + -j-F), where A and B are temperature depen- 
dent variables determined experimentally (4.097xl0 -5 
and 4.907xlCT 10 m/V for STO at 4.2K, respectively). 
This however produces a carrier density which has a sin- 




FIG. 4: Simulated 3-D carrier density as a function of 
depth for various electric potentials (OV - solid black, 100V - 
dashed red, and 200 - dotted blue) for the LAO/STO inter- 
face. This describes the increase in the overall carrier density 
at the interface as the electric potential is increased. As this 
occurs, the relative interfacial region reduces in depth. The 
inset shows the effective Fermi energy Ef as a function of 
electric potential. A fitting of this data to a quadratic poly- 
nomial provides a general functional form. 
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where Cj (i = 1, 2, and 3) are determined from experi- 
mental fits. From this, we find that 
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(8) 

where C* 4 = 4CiC 3 -C 2 Z . 

To examine the modulation of the carrier density along 
the LAO/STO interface, we compare to data provided 
by Ref. Here, a gate voltage is used to obtain 

an electric field, where F=F + V/d, V is the electric 
potential and d is the sample thickness (0.5mm). The 2- 
D carrier density is also related to the capacitance C(V) 
by 



C(V) = eS 



( d8n 



2<i 



\ dV 



(9) 



where S is the area of capacitor and 5ri2d = n 2 d(V) ~ 
n 2 d(200V) (measuring V in Volts). The cross-sectional 
area of the capacitor is determined by adjusting the sim- 
ulated peak capacitance to match the the observed ca- 
pacitance. Using this method, the cross-sectional area is 
determined to be about 5 mm 2 . Therefore, the capaci- 
tance will be 



C(V) 



Ad 



C 3 #r Fn + K 



v Ci + C 2 a(Fq + -j) 1 - 3 va. 

(10) 

Comparing to the data in Ref [15j | and adjusting only 
the Ci parameters, Figure [3] shows the capacitance and 
2D carrier density as a function of gate voltage for STO. 
Through fitting of the Sn 2 d data from Ref.[15[, the ex- 
tended parameters (Ci, C 2 , and C3) arc determined to 
be 4.25, -0.37, and 0.29, respectively with a value of Fq = 
1.2xf0 5 V/m. Here, we find that the change in the 2-D 
carrier density shifts from electron to hole doping around 
30 V. This produces the corresponding peak in the ca- 
pacitance. Further analysis in Section V will correlate 
the capacitance to the change in T c . 

To understand the accumulation of charge at the 
LAO/STO interface as a gate voltage is applied, we ex- 
amine the 3D carrier density as function of depth. The 
depth profile of the charge accumulation in STO is given 
by summing over the density profiles |£i,j(z)| 2 of all oc- 
cupied states 



n 3d (F,z) = ^2 
j=l,h 
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2ttH 2 



^(A J B lJ )e(A J B iJ )|^ J (2)| 5 



(11) 
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where AEij = Ep — Eij, and the sum over j = l,h is 
over the light and heavy electronic bands in STO. The 
inner sum is over the occupied states labelled i,j for the 
energy of eigenstate Eij and Fermi energy Ep, where 
the 8 function assures that all states Eij < Ep are 
included, m* is the effective mass (m ; * = 1.2too and 
m* t = 4.8too), mo is the electron mass, and gj = 1,2 is 
the band degenerac y 16 ' 32 The confining potential is well 
represented by a triangular shape and the corresponding 
eigenstates are therefore normalized Airy functions: 



E. 



ho 
eF 



(12) 



where n c is the normalization constant in units of 
(length)" 1 , Ek, s = f (f (i-i)) 2/S and aj = 

(2m*eF/H 2 ) 1/3 . 

Figure [4] shows the 77,3d electron density as a function 
of z. Here, the the Fermi energy Ep(V) is determined 
self-consistently by calculating 



n 2d {V) = / n 3d (z)dz 



E F (V) 



N 2d (E,V)dE (13) 



where N 2d (E, V) is the density of states 



N 2d {E,V) = 



9 j 1 



2nh 2 



(14) 



and compare to the value from Eq. [8] We can also deter- 
mine the average 713d by (77.3d) = / n\ d dz/ J n 3d dz and the 
average thickness d av — J n 3d zdz/ J n 3d dz. To compare 
the DOS to the measured capacitance, we consider 



Sn 2d {V) 
SV 



N 2d (E F {V),V) 



SE F {V) 
SV 



Ef(V) 



(15) 
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Here, it is shown that as the gate potential is increased, 
the carriers in the system are increased at the interface. 
This increases the overall density of carries and shortens 
the interfacial depth. By using the above determined 
Fermi energy Ep(V), the density of states in 2d may be 
related to the derivative of n 2d (V) with respect to V and 
hence to the capacitance as 



N 2d (E F (V),V) 



dn 2d f E ^ S 



dV 



SV 
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^S C ( V )-J ^N 2d (E,V)dE 



dE F (V) (is) 
dV 



dE F {V) 



The inset of Fig. @] shows Ep as a function of V. Us- 
ing a quadratic polynomial fit to the data, we find that 
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FIG. 5: Calculated superconducting transition temperature 
T c as a function of electric potential for the LAO/STO in- 
terface. Here, it is shown that due to the increase in carrier 
density, there is a distinct change in T c . Eventually, at some 
critical potential, the Coulomb repulsion takes over and T c is 
reduced. Inset shows the experimentally observed change in 
T c from Ref. [3] (dashed blue) and Ref. pH (dotted red). 
The difference in peak width and level can vary depending on 
the pair interaction's dependence on V as well as the exact 
determination of 8Ef(V)/ dV. 



E F {meV) = 4.30+ 0.0410 V-5.05xlCr 5 y 2 . Therefore, we 
can used the derivative to find an approximation to dV 
j dEp(V). The integral over N 2d (E, V) is dependent on 
the electric-field dependence of the effective masses and 



E, 



However, it is assumed that this contribution is 



small within the confines of the interface. 



V. EFFECT ON SUPERCONDUCTING 
TRANSITION TEMPERATURE T c 

To explore the consequences of the change in carrier 
density induced by an applied electric field on the su- 
perconducting transition temperature. We exploit the 
fact that the change in the DOS will modify the (di- 
mcnsionless) pair coupling. We conjecture that the pair- 
ing is similar to the one observed in bulk doped STOJi 
There weak coupling BCS theory provides a quantitative 
description of the observed transition temperature as a 
function of doping. Within weak coupling theory, the 
critical temperature T c below which a material becomes 
superconducting is given by 



(17) 



where Tp> is the Debye temperature (513 K for LAO/STO 
multilayers), X p is the effective interaction , and N(Ep) 
is the electronic density of states at the Fermi level in 
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d = 3 dimensions. In the case of bulk doped STO, a de- 
tailed theoretical study of the transition temperature in 
the framework of Eliashberg's theory has been performed 
Ref. [33[. The conclusion was that inter- and intra- 
valley electron-phonon processes as well as Coulomb re- 
pulsion contribute to determining T c . These calculations 
accounted quantitatively for the observed maximum of T c 
as a function of carrier density. The initial increase of T c 
at low density is connected with the increase in the den- 
sity of states, while the decrease beyond the maximum 
may be traced to the weakening of the attraction caused 
by enhanced screening. Comparing the values of density 
77.3d found above with the density at the maximum of 
T c in Ref. (33 



10 cm -3 we conclude that our 
systems are on the low density side, where density de- 
pendent screening effects of the pair interaction are less 
relevant. We therefore assume that the pair interaction 
X p may be taken as a constant, which we deduce from 
the dimensionless coupling, \ p N(Ep) = 0.13 and the 
measured DOS N(E F ) = 2 x 10 47 J -1 m _3 at density n™ 
as X p w 6.5 x 10~ 49 Jm 3 . However, the pair interaction 
may also have a dependence on the applied potential. As 
shown above the 2d electronic density of states can be di- 
rectly measured through the capacitance, and the depen- 
dence of the Fermi energy on applied voltage, while the 
latter can be inferred from the observed density-voltage 
relation. Therefore, the 3d DOS can be approximately 
expressed as 



N(E F {V),V) 



1 



Sed av j 



C(V) 



dV 



dE F (V) 



(18) 



where d av is the average thickness of the electron layer. 
There is a second term that we discuss in the supplemen- 
tary material, but as discussed above, this term is small 
and can be ignored. Hence, a change in the capacitance 
with applied electric field provides a direct effect on T c 
within this system. Figure [5] shows the calculated T c as a 
function of V for the LAO/STO heterostructure, where 
7 w 0.0045 has been adjusted to place the peak value of 
T c at the experimentally determined 0.37 K. This details 
the subtle significance of a gate voltage on the super- 
conducting state. Inset of Fig. [5] details and compares 
the experimentally observed change in T c from Ref. (l6j 
(dashed blue) and Ref. [l5[ (dotted red). 

The observed transition into the superconductive 
state is most likely of the Kosterlitz-Thouless-Berezinski 
(KTB) type. In general the transition temperature Tktb 
is related to the bulk transition temperature T c by 3 -! 



Tktb 
T, , 
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KTB 



2.18 



Rc 



R 



sheet 



(19) 



where R s heet is the sheet resistance and Rq = h/e 2 m 
4.12fcf2 is the resistance quantum. Here the function / is 
defined in terms of the gap parameter A(T) as 



f 



A(T) 
A(0) 



tanh 



A(T) 
IT 



(20) 



In the limit of weak disorder, Tktb ~ T c . At strong 
disorder, when R s heet >> 2i?g one gets a suppres- 
sion of the KTB transition temperature as Tktb = 
2.18T C R Q /R 

sheet ■ 



VI. CONCLUSION 

In conclusion, we examined the effect of an electric 
potential on the carrier densities and superconducting 
transition temperature for the conducting LAO/STO in- 
terfaces. We find that the carrier density and supercon- 
ducting temperature can be significantly modified by an 
external electric field. Recently, there have been experi- 
mental observations that have explored the relationship 
between an applied potential on STO heterostructures 
and the increased carrier density or capacitance. Here, 
we take this relationship further and assume it is the field 
effect that, through nonlinear capacitance, changes car- 
rier density. That in turn affects the superconducting 
coupling and ultimately T c . We present a mechanism of 
this phenomenon assuming a conventional pairing state 
at the interface, which shows a possible route to tunable 
material properties and suggests an exciting perspective 
on the way to tunable superconductivity. 
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